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Abstract 



The extension of the first-principles generalized pseudopotential theory 

•rH . 

(GPT) to transition-metal (TM) aluminides produces pair and many-body 
interactions that allow efficient calculations of total energies. In aluminum- 
rich systems treated at the pair-potential level, one practical limitation is a 
transition-metal over-binding that creates an unrealistic TM-TM attraction 
at short separations in the absence of balancing many-body contributions. 
Even with this limitation, the GPT pair potentials have been used effectively 
in total-energy calculations for Al-TM systems with TM atoms at separations 
greater than 4 A. An additional potential term may be added for systems 
with shorter TM atom separations, formally folding repulsive contributions 
of the three- and higher-body interactions into the pair potentials, result- 



ing in structure-dependent TM-TM potentials. Towards this end, we have 
performed numerical ab-initio total-energy calculations using VASP (Vienna 
Ab Initio Simulation Package) for an Al-Co-Ni compound in a particular qua- 
sicrystalline approximant structure. The results allow us to fit a short-ranged, 
many-body correction of the form a(ro/r) b to the GPT pair potentials for Co- 
Co, Co-Ni, and Ni-Ni interactions. 
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I. INTRODUCTION 



Total-energy calculations are an important tool in theoretical condensed-matter physics, 
giving insight into structures and mechanical properties of solids [0,0. Accurate calcula- 
tions of total energy are notoriously difficult, however. Theoretically, one must solve the 
Schrodinger equation simultaneously for all electrons in the presence of fixed atomic nu- 
clei. Density-functional theory (DFT) simplifies this problem by reducing it to the 
self-consistent solution of Schrodinger's equation for a single electron in a potential that de- 
pends upon the electron density. Even with this simplification, full ab-initio DFT electronic- 
structure methods are computationally demanding ||, usually limited to systems of less than 
a hundred atoms, and may not yield immediate physical insight once an answer is obtained. 

Instead, one may expand the total energy in terms of pair and many-body interatomic 
potentials |l||JJ§-P], so that the total energy appears as an explicit function of atomic sepa- 
rations. Depending on the physical system under study and the type of information sought, 
the expansion may often be truncated after a small number of terms. Such a truncated ex- 
pansion trades off a degree of accuracy in favor of computational simplicity and potentially 
greater physical insight as compared with a full ab-initio electronic-structure approach. 

Many metallic systems been studied using such quantum-based interatomic potentials, 
including aluminum and its alloys with both transition and non-transition metals PHHfl- 
These potentials are especially simple in the case of non-transition metals. There, the d- 
electron energy bands are either empty or else are completely filled and deeply buried below 
the Fermi energy level, allowing rapidly convergent expansions of the total energy and an 
accurate description in terms of only radial-force interactions . The presence of partially 
filled ci-bands at or near the Fermi level in transition metals complicates the analysis. The 
occupied ci-band electronic states are highly localized in the vicinity of the atoms, leading 
to directional or covalent bonding with a strong angular dependence. Consequently, total- 
energy expansions will not converge as quickly as for non-transitions metals, and three- and 



higher-body angular-force interactions may contribute significantly Bfl|l6|.|T7|j 
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Moriarty [[17] has developed a rigorous DFT treatment of interatomic potentials for 
transition metals (TMs) in the context of the generalized pseudopotential theory (GPT). 
The treatment was later extended to binary and ternary alloys of aluminum with first row 



transition metals |H |. These studies found that three- and four-body interactions could be 
important in determining thermodynamic and mechanical stability of structures with large 
TM concentrations. The explicit treatment of <i-electron interactions in the GPT produces a 
strong attractive interaction at unphysically short distances in the pair potentials, which is 
balanced by repulsive forces contained in the many-body interactions. For specific structural 
environments, however, it is possible to directly modify the short-ranged part of the TM- 
TM pair potentials to remove this unphysical attraction, so that a truncation of the total 
energy expansion at the level of pair potentials will be more accurate when transition metal 
atoms are near neighbors. Here we wish to consider the construction of such effective pair 
potentials for important transition-metal aluminide systems. 

One motivation for this study is the need for fast total-energy calculations in Al-TM 
systems with short TM-TM separations to enable structural relaxation, and more generally, 
molecular-dynamics and Monte-Carlo simulations. We focus our attention on Al-Co-Ni com- 



pounds in decagonal quasicrystalline structures flT8[] . The precise modifications required in 
the Co-Co, Co-Ni, and Ni-Ni pair potentials depend on the particular structure studied, but 
they should be at least approximately valid for many similar structures. Furthermore, the 



modifications obtained may allow us to treat Al-Co-Cu and Al-Cu-Ni decagonal phases [19 



because the Cu-Cu interactions do not appear to require modification ||14|| . Limited numbers 
of ab-initio electronic-structure calculations, which effectively sum the pair and many-body 
total-energy contributions, are sufficient to determine the required modifications, and this 
is the strategy that we follow here. 

We intend to apply these potentials to predict the structures of decagonal quasicrys- 



tals |2(J . A great deal of experimental data is available that identifies the positions of most 
atoms and identifies the chemical identity of many of those. However, in order to determine 
the quasicrystal structures from X-ray diffraction, one faces degenerate structures because 



elements near each other in a row of the periodic table (such as Co, Ni and Cu) have similar 
X-ray form factors. A common approach to this problem is to supplement the experimental 
data with total-energy calculations. This approach is well established in crystallography [fZL] . 

The effective pair potentials developed here can be applied to total energy calculations in 
quasicrystals and related structures with a great reduction in computational times compared 
with the ab initio electronic-structure calculations themselves. The time savings results from 
two general features of the potentials. First, for a given atomic volume and composition 
the potentials may be precalculated and then applied repeatedly with a simple lookup and 
interpolation. Second, to calculate the change in energy when a single atom is moved, only 
interactions affecting that atom are needed. If the interactions are cut off at a finite spatial 
separation, the time required to calculate the change in total energy becomes independent 
of the number of atoms in the complete structure. This is so-called order-A scaling. In 
contrast, ab-initio electronic-structure methods must recalculate the entire system when a 
single atom is moved, typically resulting in order- A 3 scaling. 

In Sec. |TJ, GPT interatomic potentials are briefly reviewed and the issues surrounding 
the truncation of the total energy expansion at the pair-potential level in Al-TM systems are 
discussed. Section |TJ gives details about the scheme we employ to determine the needed 



modifications to the TM-TM pair potentials. In Sec. [TV], we discuss the results of our full ab- 
initio electronic-structure calculations and the modified TM-TM pair potentials developed 
using them. 



II. GPT INTERATOMIC POTENTIALS 

The generalized pseudopotential theory starts with a full ab-initio DFT representation 
of the total energy in the standard local-density approximation (LDA). The usual small- 
core approximation is used to separate the treatment of valence and core electrons and the 
electron-ion interaction for the valence electrons is handled by means of optimized nonlocal 
pseudopotentials. A mixed valence-wavefunction basis is employed allowing sp states to be 
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represented as superpositions of plane waves, while d states are represented in terms of local- 
ized, atomic-like d states. The electron density and total energy are systematically expanded 
in terms of the resulting weak matrix elements in this basis: sp pseudopotential matrix el- 
ements Wkk'i sp-d hybridization matrix elements Akd and Skd, and d-d tight-binding-like 
matrix elements A^' and S^d 1 - I n real space, the total energy may be cast in terms of 
well-defined interatomic potentials, which can be calculated as functionals of these matrix 
elements. For a general multicomponent alloy, the GPT total-energy expansion takes the 
form 



E tot (Ri) = NE vol (n, c) + i E ' vf(R& ^ c ) + 1 E E ' Rjk, R ki ] n, c) + 

a(3 ij a/3"/ ijk 



(1) 



where R is the set of all positions of N ions in the metal, E vo \ is a volume term that includes 
all collective and one-ion contributions that are independent of structure, and , 1 , . . . 
are the two-, three-, and higher many-ion interatomic potentials. The primes on sums over 
ion positions exclude all self-interaction terms. Indices cn,/3, 7, ... run over all chemical 
species, and indices k, . . . run over the individual ion positions. The volume term and 
all of the interatomic potentials depend on the atomic volume Q and a composition vector 
c, but are independent of structure. The potentials are functions of the relative positions 
of small subsets of atoms, independent of the positions of all other atoms in the system. 
The entire dependence on the structure comes analytically through the summations over 
ion positions. This makes these potentials transferable among different structures at fixed 
atomic volume and composition. The full details of the first-principles GPT for transition- 



metal systems are given in Refs. [O] and fiM. A simplified model version of the theory has 



also been developed [^] , using canonical <i-bands to obtain analytic representations of the 
multi-ion potentials. 

In general, the separation of the total energy into two- and higher-body terms is not 
entirely unique, since one can always add contributions to the pair potential provided 
one makes suitable subtractions from v^ 1 and/or higher-body potentials. Within the GPT, 



the uniqueness of the potentials is established by ensuring that their desired properties of 
structure independence and full transferability are consistent with the matrix elements that 
define them. In this regard, the total energy is normally calculated to second order in the 
pseudopotential Wkk', so that sp contributions enter only in the volume term and the pair 
potentials. The TM d-d and sp-d contributions to each potential, on the other hand, are 
carried to all orders in the matrix elements A^', etc. Terms are allocated to pair- and 
many-body potentials according to how many distinct ionic positions explicitly enter. Thus, 
for example, the TM pair potentials v 2 _ contain contributions that are even powers 
of Add' associated with repeated hopping of d electrons back and forth between a pair of 
ions, with the leading term proportional to A dd /A d / d . The three-ion TM potentials contain 
corresponding terms of third order proportional to A dd > A d ' d »A d n d and terms of fourth order 
proportional to A 2 dd ,A 2 d , d „, as well as higher-order terms. The four-ion TM potentials start 
at fourth order in A dd i. 

The tight-binding-like d-d contributions to the TM potentials are modulated by sp-d hy- 
bridization, d-state nonorthogonality, and other factors such as <i-band filling, but nonethe- 
less, they give valuable insight into the expected short-range behavior of the potentials. In 

rp TlT rp i\/r 

particular, one expects a strong attractive contribution to v 2 ~ at short distances result- 
ing from the second-order term in A dd > . This term is attractive because it directly relates to 
the second moment of the ci-band density of states and hence to the <i-band width and the 
additional cohesion provided by partial rf-band filling. The attraction is strong at short dis- 
tances because the matrix element for atoms separated by distance r varies roughly as r~ 5 , 
the behavior obtained for pure canonical d bands. In addition, one expects this attractive 
contribution to be maximum near half-filling of the d bands and to vanish for completely 
filled or empty d bands. Thus the expected overbinding in v ™-™ show a clear chemical 
dependence, with decreasing magnitudes for Co-Co, Co-Ni, and Ni-Ni interactions. 

For short-range TM interactions, repulsive contributions from higher-order terms in Add' 
will balance the attractive contribution of the second-order term in t>J M ~™, provided that 
the local concentration of TM atoms is sufficiently high. In general, the detailed balance 
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obtained involves the three-, four-, and possibly higher-ion potentials. Near half-filling of the 
d bands, however, the repulsive contributions will be dominated by the fourth-order terms 
in v 3 ~ ~ . This in turn suggests a simple scheme to modify the TM pair potentials 
at short range to incorporate the balance directly, a scheme that we will develop in Sec. 
First, however, we examine the actual calculated GPT pair potentials for the Al-Co-Ni 
system of interest here. 

Figure [l] shows the Al-Al and Al-TM pair potentials for Al-Co-Ni [fL4| . These are calcu- 



lated in the aluminum- rich limit, but in practice they do not depend strongly on composition. 
The first minima of the Al-TM pair potentials occur near 2.3 A with depths of about 0.2 eV 
(Al-Ni) and 0.3 eV (Al-Co). Rather than a potential minimum, the Al-Al potential exhibits 
a shoulder near 3 A. The TM-TM pair potentials are shown in Fig. 0. As expected, the 
TM overbinding is most evident for Co. The first minimum in the Co-Co potential has a 
depth of 2.1 eV at 1.7 A. The corresponding Ni-Ni potential depth of 0.1 eV at 2.2 A is 
not obviously problematic, but in the following we will find it requires some modification. 
In the present applications, the Co-Ni pair potential V2° Nl is defined as an average of the 
Co-Co and Ni-Ni potentials, 



CoNi _ („,CoCo , NiNi 



+ <^)/2. (2) 

This amounts to a perturbative expansion of in the difference in atomic number Z a — Z^. 
Clearly, V2° N% so-defined suffers overbinding due to the overbinding of v^ 000 ■ 

We wish to devise effective pair potentials for Al-Co-Cu and Al-Cu-Ni as well as Al- 
Co-Ni. Previously, the Al-Cu potentials were found to be well behaved up to large Cu 
composition ||14j| , so no modification of is suggested. Our modification to V2° Co ob- 

tained for Al-Co-Ni compounds may be approximately valid for these other compounds. We 
previously defined |H| V2° Cu as equal to v% because Ni lies between Co and Cu in the 



periodic table. Thus our modified Ni-Ni potential should serve as an approximate modified 
Co-Cu potential. For the modified Cu-Ni potential we may take (v2 UCu + v^ Nl )/2, using 
the modified Ni-Ni potential. The Al-Co-Cu and Al-Co-Ni potentials so-obtained will, of 



course, still need to be validated using full ab-initio calculations. 

III. MODIFICATION OF PAIR POTENTIALS 

As discussed above, the short-ranged attraction in the TM-TM pair potentials is bal- 
anced by repulsive terms contained in the three- and higher-body potentials. If one chooses 
to truncate the GPT expansion at the pair potential level, these repulsive many-body contri- 
butions must be "folded" into effective pair potentials. Formally, one may define an effective 
pair potential by averaging over atomic positions, holding a single pair of ions fixed: 

v *ff = v f + < v «fh > + < >+.... (3) 

Such potentials have been previously considered in the context of the simplified model 
GPT []22] and canonical d bands for central transition metals. There it was found that 
the four-body interaction oscillates with respect to angles between atoms, with a nearly 
zero mean, so it does not contribute significantly to v% . The third-order contribution to 
v% also approximately averages away, but the fourth-order contribution to v 3 ^ 7 contributes 
strongly, yielding a short-ranged repulsive term proportional to A 2 dd ,A 2 d , d ~ r~ 20 balancing 
against the attractive second-order term in v% . 

Inspired by the short-ranged repulsion found in Eq. (§) and the power law variation of 
Add' within the model GPT ||22|| , we propose to modify the full GPT pair potentials v% by 
adding terms of the form 

U a0 (r) = a(r /r) b , (4) 

where a and b are positive and depend upon the elements a and (3 of the pair potential 
modified. Our expectation, which is confirmed below, is that b is large in all cases, so that 
U a/3 is indeed short-ranged. In our applications, the quantity tq is taken as a common atomic 
separation in quasicrystals of 2.55 A. Then at a fixed atomic volume and composition the 
effective pair potential can be written as 
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V a/3 (r) = vf(r) + U a ?(r) . (5) 

We determine the unknowns a and b by matching energies and forces obtained from full 
ab initio electronic-structure calculations on a quasicrystal approximant. Cockayne and 



Widom p3||24| previously suggested a structure for decagonal Al-Co-Cu. An approximant 
of that structure is shown in Fig. |3] with Ni atoms replacing Cu. The orthorhombic unit 
cell (a=23.3 A, b=7.57 A, c=4.09 A) contains 50 atoms ( A^CoioNie) . Most atoms occupy 
either z=0.25 or z=0.75 layers. Al atoms at the centers of hexagons occupy the z=0.5 layer. 
Two Co atoms occupy symmetric positions around these central Al atoms. In Al-Co-Cu, 



alternation of Co and Cu on tile edges is thought to be energetically advantageous [ 2J] . We 
find that alternation of Co and Ni shown in Fig. [3] is slightly disadvantageous in Al-Co-Ni. 

To investigate TM bonding energetics, we alter the basic structure shown in Fig. ^ by 
swapping a Co atom on a horizontal tile edge (atom b in Fig. |||) with the Ni atom on the 
other horizontal tile edge (atom c). Focusing on near-neighbor interactions, we find this 
swap of atoms replaces four Co-Ni bonds with two Co-Co and two Ni-Ni bonds, all of length 
= 2.55 A. These numbers are twice as large as is apparent by inspection of Fig. |3|, the extra 
factor of two coming from periodic boundary conditions in the direction perpendicular to 
the plane. 

Now consider the energy change evaluated using pair potentials. Atoms b and c occupy 
nearly equivalent sites. An exact symmetry in the Al atom positions guarantees that no 
bond involving an Al atom is affected by the swap. We already noted the change in TM 
near-neighbor interactions. At further neighbors, with separations of 4.6 A or greater, we 
also find interchanges between Co-Co and Ni-Ni bonds for pairs of Co-Ni bonds. If the 
approximate form (|2|) of v 2° Nl as the average of v^ oCo and were valid, all changes in 

bonding would exactly cancel each other, resulting in a vanishing energy change. We presume 
that approximation (|2|) is more accurate at large separations than small separations. Thus 
we attribute the entire energy change of the be swap to near-neighbor energy differences 

AEi = 2V CoCo + 2V Nim - 4V CoNi , (6) 
10 



where V alS denotes the strength of the pair potential evaluated at the near-neighbor distance 
2.55 A. 

Next we swap one of the Co atoms inside the tiles (atom e) with one of the Ni on 
a horizontal tile edge (atom a). Two Co-Ni bonds are broken and two Co-Co bonds are 
produced after this swap. All other interactions that are affected are Al-TM interactions, 
which we presume to be described accurately by the GPT pair potentials. This swap energy 
can be written as: 

AE 2 = 2V CoCo - 2V CoNl + V Al ™ , (7) 

where V Al ™ represents a calculable collection of interactions between Al atoms and TM 
atoms at many separations. V Al ™ should be described accurately by the unmodified GPT 
pair potentials. 

Lastly, we replace the Co-Ni pair on one horizontal tile edge (atoms c and d) with Al 
atoms. Then we swap one of the newly introduced Al (at position c) with a Ni atom on the 
other horizontal tile edge (atom a). This breaks two Co-Ni bonds. All other interactions 
are either Al-TM or Al-Al interactions, and again those are described well within the GPT. 
The energy change of this swap is 

AE 3 = -2V CoNi + V Al ™ + V AlAl , (8) 

where V AlAl and y Al ™ represent collections of interactions involving Al atoms that, as 
before, we presume to be accurately calculable within the unmodified GPT. 

Full ab initio values for the energy changes AEi, AE 2 and AE 3 were calculated using 
VASP (Vienna Ab Initio Simulation Package) p5[ . VASP calculates total energies within 
the local-density approximation using pseudopotentials to treat valence-core electron in- 
teractions. We performed calculations using a 4x4x4 k-space grid and also using a 4x4x8 
k-space grid to observe the convergence as k-points are added. All calculations were done 
using medium precision which is expected to be sufficient for our needs. We iterate the 
self-consistent calculation until an accuracy of 10~ 6 eV is achieved. 
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By comparing the energy differences AEi, AE2 and AE3 calculated by VASP with 
the same quantities calculated with the unmodified GPT potentials, we can obtain the 
values of U af3 evaluated at the near-neighbor separation 2.55 A. Specifically, when energy 
changes calculated by the unmodified GPT are subtracted from energy changes calculated 
by VASP, assuming that the contributions V AlAl and y M ™ are accurately calculated with 
the unmodified GPT, we find 

AE\ ASP - AE? PT = 2U CoCo + 2U mm - AU c ° m 

AE^ ASP - AE% PT = 2U CoCo - 2U mm (9) 
ae vasp _ ae gpt = _ 2U Com 

Since each correction U al3 (r) involves two unknowns, a and b, Eq. (|9|) consists of three 
equations in six unknowns. Additional information is obtained from the forces on atoms 
calculated by VASP. By examining the forces on the Co-Ni pair (atoms c and d) in Fig. |3|, 
and on the Co-Co and Ni-Ni pairs created by the be swap, we obtain three additional 
equations governing the derivatives of U a ^ at the near-neighbor separation. This additional 
information allows closure of the equations and determination of the unknowns. 



IV. RESULTS 

Table [I] shows the energy differences AEi in Eqs. ©-© calculated using GPT pair 
potentials and VASP. Comparing the VASP data for the two grid sizes, we note that the 
signs and approximate magnitudes of AEi are consistent with each other. One immediate 
result from Table | is that mixed Co-Ni bonds are disfavored over pure Co-Co and Ni-Ni 
bonds. The energy difference AE\ results from breaking four Co-Ni bonds and producing 
two Co- Co and two Ni-Ni bonds. AE\ calculated by VASP is negative, showing that the 
swap lowers the system energy. This means that for Al-Co-Ni, similar TM atoms prefer to 
reside near each other on the tile edges. Cockayne and Widom found the opposite for the 



case of Al-Cu-Co using mock ternary potentials |2q], and this was confirmed later using a 
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full ab-initio technique p4 



Also concerning the calculated values of AEi, we see that the averaged potential approx- 
imation (0) is fairly accurate. GPT yields AE\ = because it employs this approximation 
here. The small value of AEi obtained by VASP confirms that this approximation is not 
far off the mark. 

Figure § shows the x-component of the total force on certain TM atoms. Our (4x4x4) 
and (4x4x8) VASP calculations yield forces that agree to 0.06 eV/A or better. We examine 
the horizontal bonds ab and cd in Fig. ^| in both the original and swapped configurations. 
As expected, at 2.55 A, GPT pair potentials predict attractive forces between TM pairs 
while the actual forces obtained from VASP are repulsive. The small force asymmetry on 
atoms in the Co-Ni pair is due to the different ways Co and Ni atoms interact with their 
surrounding environments. The difference between the forces calculated by our two methods 
is greatest for Co-Co bonds and smallest for Ni-Ni bonds, consistent with our expectation 
that overbinding is more severe for Co than for Ni. 

Calculated modifications to the GPT pair potentials are given in Table [TI]. Examining the 
magnitude of U af3 at r=r =2.55 A (i.e., the value of a), we note that JJ NtNl is smaller than 

jjCoCo 

as is expected since Ni is closer to the end of the 3c? transition series, with its <i-bands 
almost full. It should be noted that r=2.55 A is not the potential minimum. It is rather 
the nearest-neighbor distance at which the calculations were performed. The quantities V af3 
and F af3 are, respectively, the energy and force calculated from the modified GPT potentials 
[Eq.( ^|)] at the near-neighbor distance r . The large values of b we obtain show that our 
modifications of the GPT pair potentials fall off rapidly beyond the near-neighbor separation 
and confirm our expectations based on Eq. (3). The modified potentials are illustrated in 
Fig. [|. The (4x4x4) and (4x4x8) VASP calculations agree in positions of the potential 
minima to about 0.05 A and agree in the values at the minima to about 0.02 eV. 
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V. DISCUSSION 



The original GPT interatomic potentials were derived from first principles without refer- 
ence to specific structures. Their applicability to, and transferability among, a broad range 



of structures was verified [13]. The pair potentials alone apply to Al-rich structures in which 
TM atoms are well separated, but fail due to an unphysical short-ranged TM attraction. 
Our modification of the GPT potentials is restricted to TM pair potentials. Since the correc- 
tion U af3 (r) is negligible beyond 3 A, our corrections only affect energies and forces among 
neighboring TM atoms. 

In idealized decagonal AINiCo and AlCuCo quasicrystal models, neighboring TM atoms 
always occur in specific atomic environments consisting of zig-zag chains of TM atoms at 



2.5 A spacing surrounded by Al atoms at special positions | 20y23|j24 |. A variety of qua- 
sicrystal and approximant structures differ in the arrangement of these chains in space, but 
share the same local structure around the TM atoms. Our modified GPT potentials are 
therefore strictly transferable within this class of quasicrystal structures. 

While the modified GPT potentials are not strictly transferable outside this special 
class of structures, we do believe they are qualitatively transferable. Comparing figures 
2 and 5, we have replaced obviously unphysical pair potentials with a set that appears 
qualitatively realistic. Both original and modified pair potentials are transferable among 
structures without TM neighbors. The modified pair potentials provide a reasonable, though 
non-rigorous, extension to structures with TM neighbors. 

As a test of our modified pair potentials, we relaxed the basic structure using both 
original and modified GPT pair potentials, and using VASP. In all three runs we relaxed 
the structure until all atomic forces were less than 0.001 eV/A. Under VASP relaxation (4x 
4x4 k-space grid), Al atoms moved 0.19 A on average, followed by Ni with an average 
displacement of 0.16 A and Co with an average displacement of 0.15 A. Relaxation under 
the modified GPT potentials produced a structure close to the VASP relaxed structure. 
The differences between modified GPT and VASP relaxed positions are less than 0.07 A for 
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every TM atom, with an average difference of 0.05 A/TM atom. In contrast, relaxation 
under the original GPT pair potentials yielded a maximum TM relaxed position difference 
of 0.20 A and an average TM difference of 0.13 A. 

As a result of our modifications, the TM atoms relax in the correct directions and move 
approximately the correct distances. This is not the case using the original GPT potentials. 
Since TM separations move away from r under relaxation, our proper behavior under 
relaxation demonstrates a type of transferability of the potentials. 
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TABLES 

TABLE I. Total-energy differences denned by Eqs. ©-(H) as calculated by VASP and GPT. 
Units are eV/cell. 



Energy 


GPT 


VASP (4x4x4) 


VASP (4x4x8) 


AEi 


0.000 


-0.020 


-0.031 


AE 2 


0.116 


0.298 


0.279 


AE 3 


-0.945 


-1.384 


-1.419 



TABLE II. Modifications for GPT transition-metal pair potentials, U a/3 (r) = a(ro/r) b where 
ro=2.55 A. The quantities V al3 and F a/3 are the energy and force calculated at r=ro from the 
modified GPT potential, Eq. (Hj). Units of a and V a @ are eV while b is dimensionless and F a @ has 
units of eV/A. 



a(3 


a 


b 


yap 


pot/3 


CoCo 


0.319 


16.6 


0.0946 


0.978 


CoNi 


0.237 


19.3 


0.0941 


0.994 


NiNi 


0.140 


21.3 


0.0779 


0.674 
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FIGURES 




FIG. 1. GPT interatomic pair potentials for Al-Al, Al-Co and Al-Ni. 
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FIG. 2. GPT interatomic pair potentials for Co-Co, Co-Ni and Ni-Ni. 



20 




21 



VASP GPT 

Co Co Co Co 




Ni Ni Ni Ni 




Ni Co Ni Co 




FIG. 4. Horizontal components of forces (in eV/A) on transition metal atom pairs calculated 
from the GPT and VASP. 
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FIG. 5. Modified transition-metal GPT pair potentials using the parameters in Table ||. 
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